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We propose a formulation of an absorbing boundary for a quantum particle. The formulation is based on a Feynman- 
type integral over trajectories that are confined to the non-absorbing region. Trajectories that reach the absorbing 
wall are discounted from the population of the surviving trajectories with a certain weighting factor. Under the 
assumption that absorbed trajectories do not interfere with the surviving trajectories, we obtain a time dependent 
absorption law. Two examples are worked out. 
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> 

m ■ 
O ■ 

oo : 
o . 
■ 

ON 1 

9 L,: 
^— > ■ 

G ■ 

3 : 



The purpose of this letter is to propose a Feynman-type 
integral to describe absorption of particles in a surface 
bounding a domain. The need for such description arises 
for example in scattering theory, in the description of a 
photographic plate, in the double slit experiment, in neu- 
tron optics, and so on. The optical model is often 
used to describe absorption in quantum systems. This 
model is based on analogy with electro-magnetic wave 
theory. It is not obvious that the methods of describing 
absorption in Maxwell's equations carry over to quan- 
tum mechanics because the wave function of a particle 
does not interact with the medium the way an electro- 
magnetic wave does. In particular, in classical quantum 
theory, unlike in electromagnetic theory, the wave func- 
tion does not transfer energy to the medium. 

The main tenet of our Feynman-type integral descrip- 
tion of absorption is that trajectories that propagate into 
the absorbing surface for the first time are considered 
to be instantaneously absorbed and are therefore termi- 
nated at that surface. The population of the surviving 
trajectories is therefore discounted by the probability of 
the absorbed trajectories at each time step. 

The process of discounting can be explained as fol- 
lows. In general, if the trajectories are partitioned into 
two subsets, the part of the wave function obtained from 
the Feynman integral over one cannot be used to cal- 
culate the probability of this subset, due to interference 
between the wave functions of the two subsets. How- 
ever, in the physical situation under consideration, such 
a calculation may be justified. Our procedure, in effect, 
assumes a partition of all the possible trajectories at any 
given time interval [t, t+ At] into two classes. One is a 
class of bounded trajectories that have not reached the 



surface by time t + At and remain in the domain, and 
the other is a class of trajectories that hit the surface for 
the first time in the interval [t,t+ At]. We assume that 
the part of the wave function obtained from the Feynman 
integral over trajectories that hit the surface in this time 
interval no longer interferes with the part of the wave 
function obtained from the Feynman integral over the 
class of bounded trajectories in a significant way. That 
is, the interference is terminated at this point so that 
the general population of trajectories can be discounted 
by the probability of the terminated trajectories. This 
assumption makes it possible to calculate separately the 
probability of the absorbed trajectories in the time inter- 
val [t,t+ At]. 

The discounting process constitutes a coarse-graining 
procedure for a large quantum system, describing the 
absorber. The assumptions we make can be viewed as 
the mathematical expression of quantum irreversibility 
of absorption, because the absence of interference sepa- 
rates the two classes of trajectories for all times. Thus 
the trajectories that stop to interfere can be discounted 
from the population of trajecotires inside the domain. 
Our model is only one aspect of irreversible processes. 

Under the above assumptions the discounting proce- 
dure leads to a Feynman-Kac integral with a killing mea- 
sure which in turn leads to a Schrodinger equation 
with zero boundary conditions on the absorbing surface 
and complex valued energy which depends on the wave 
function of the class of bounded trajectories at each time 
t. 

We obtain a decay law that is not an exponential rate, 
in general. It reduces to a rate only for a single energy 
level initial condition. In this case, the wave function is 
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a solution of the Schrodinger equation with infinite walls 
at the absorbing boundaries and constant complex poten- 
tial that depends on the initial energy. If, however, there 
are two initial energy levels, the decay law depends on 
both energy levels and contains an oscillatory term with 
beat frequency. Similar beats in the decay law occur if 
the initial wave function contains any number, finite or 
infinite, of energy levels. 

We consider another example, of a Gaussian packet of 
free particles traveling with a given mean velocity toward 
an absorbing wall. We find that the packet is partially re- 
flected and partially absorbed and calculate the reflection 
coefficient. The reflected packet results from trajectories 
that never reached the wall. This behavior is mainly due 
to classically forbidden trajectories. In contrast to scat- 
tering, this is not the same reflection as that in a finite 
or infinite barrier, because reflection in a finite barrier 
depends on the shape of the barrier and reflection in an 
infinite barrier is totally elastic, whereas the reflection we 
obtain is discounted by a finite constant factor. 

Consider the class a a ,b of continuously differentiable 
functions x(t) for < r < oo such that a < x(t) < b for 
all < t < oo and such that x(0) = xj, x(t) = x. The 
class a a ,b consists of bounded trajectories that begin at 
xi and end at x. We define the Feynman integral over 
the class a a ^ by 



K(x,t) 



= lim a N 



where 
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exp<j -S[x(-),t]}Vx(-) 



(1) 



J ...J exp^-S(x ,...,x N ,t)j JJ dxj, 



a=l m X 1 ' 2 
I 2mhAt J 

Next, following the method of R], we show that 
K(x, t) satisfies Schrodinger's equation and determine 
the boundary conditions at the endpoints of the interval 
[a, b]. We begin with a derivation of a recursion relation 
that defines K(x,t). We set 



K N (x N ,t) 



N-l 



a N J ...J exp |^5(x , . • . ,z/v,*)j dxj, 

then, by definition, K(x,t) = limjv— »<» Kn{x, t). We 
have therefore the recursion relation 



K N (x,t) = a I exp<{ - 



m(x — xjv-i) 2 



2At 

x isTjv_i(xjv_i,tjv_i)da;jv_i. 



V(x)At 



The following derivation is formal, a strict derivation can 
be constructed along the lines of (7| . We expand the func- 
tion Kisf-x(xjsr-i,tN-i) in (||) in Taylor's series about x 
to obtain 



K N (x,t) 



(3) 



-iV(x)At/h 



exp 



K N -i(x,t N -i) - (x - x N -i) 



dK N -i(x,tN-i) 



dx 



-{x-Xn^) 2 — 

O ((x - x N ^if^j 



d 2 K N _ 1 (x 1 t N ^i) 



dx 2 
dxjf-t. 
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We evaluate the integrals in eq. (||) separately for x in- 
side the interval [a, b] and on its boundaries. This leads 
to the Schrodinger equation 



if) 



dK(x,t) 

at 



% 2 d 2 K(x,t) 
2m dx 2 



+ V(x)K(x,t) 



(4) 



for a < x < b with boundary and initial conditions 

K[a, t) = K(b, t)=0 for t > (5) 

K(x, 0) = S(x - xi) for a < x < b. (6) 

Obviously, eqs.(Q)-(||) are identical to those of a particle 
bounded by infinite potential walls. 

The same result was obtained in || for the case V(x) = 
by a different method. Our method of calculation is 
essential for calculating the Feynman integral with ab- 
sorbing boundaries. 

First, we calculate the discretized Feynman integral 
to survive (not to be absorbed) the time interval [0, At] 
and to find a trajectory in time At at a point x in the 
interval [a, b]. According to the above assumptions, the 
discretized Feynman integral for trajectories initially in- 
side the interval [a, b] that propagate to the endpoint a 
for the first time in the time interval [0, At] is 



(2) 



-0i (a, At) = a J * (a;o)exp|-S'(xo,a, At) j dx . 

Therefore, the probability density of finding a trajectory 
at the point a in the time interval [0, At] is 

l^i(a,Ai)| 2 , 

and there is an analogous expression for the probabil- 
ity density of finding a trajectory at the point b in the 
time interval [0, At] . It follows that the probability of a 
trajectory to be absorbed in the time interval [0, At] is 

P 1 {At) = X a \Ma^t)\ 2 +A b |Vi(6,At)| 2 , 

where A Q and Ah are characteristic lengths (see discussion 
at the end of the letter). Thus the discretized Feynman 
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integral to survive the time interval [0, At] and find a 
trajectory in time At at a point x in the interval [a, b] is 

*i(x,Ai) = v/l-Pi(At)a / ^o(xo) x 

J a 



exp ^ — S l (xo, x, At) j dxQ 
y/l-P^At) Kx(x,At). 



(7) 



Next, we calculate the discretized Feynman integral to 
survive the time interval [At, 2 At] and find a trajectory 
in time 2At at a point x in the interval [a, b]. Accord- 
ing to eq.(^|), given that a trajectory survived to time 
At, its discretized wave function is Ki(x, At) so that the 
discretized Feynman integral to propagate' to the point a 
is 



■0 2 (a,2At) = a J Ki(x, At) exp ^-S(x,a, At) j dx. 

Proceeding this way, we find that the discretized Feyn- 
man integral to survive the time interval [0, NAt] and 
find a trajectory in time NAt at a point x in the interval 
[a, b] is 



H> N (x,NAt) 



N-l 



x II 0--Pi(j*t))K N (x,t), 



(8) 



*0, t) = y/1 -P(t)K(x, t), (13) 
and 1 - P(t) is given by @. 
Examples 

First, we consider a particle with two absorbing walls 
at x = ±a and zero potential. We assume that A_ a = A a . 
The wave function is given by 



K(x, t) = ^ A n exp (~2^-t) sm _* 

n=l v ■* 



so that 



-K(±a,t) 



oo oo 



A k A n 2knm /^fc+n 

fc2 _ n 2 j^ fl 2 



n=l fc^n 



1 — exp 



> 2 



2ma 2 



,2 " 2 7T 2 



For a particle with a single energy level the wave func- 
tion decays at an exponential rate proportional to the 
energy. However, if there are more than just one level, 
the exponent contains beats. For example, for a two level 
system with real coefficients, we obtain 



where 

Pj(jAt) = X a \^(a,jAt)\ 2 + X b |^(6, At)| 2 . 

It remains to calculate the survival probability 

N-l 



(9) 



1 - P(t) = exp 



(A 2 k 2 + A 2 n n 2 )t- 



4mA k A n h (k 2 

~^rr Sin ■ 



h{k 



2ma 2 



1 - P(t) = j™^ IJ ^ ~~ P J ^ A *)) • ( 10 ) The strongest beats occur for k = 2, n = 1 with fre- 

quency u h2 = f^r- Setting Ai = A 2 = y/T/2 and 



It can be shown Q that the probability Pj(jAt) is given 
by 



PjijAt) 



fiAt 



2irm 



d_ 

dx 



K s -! (a, (j - l)At) 



Tx Ki-i{b, (j-l)Ai) 



so that eq.(|Io|) gives 
l-P(t) = 



exp ■ 



A, 



dx 



K(b, t) 



dt. 



(11) 



(12) 



Now, it follows from eqs.(||) that the wave function of the 
surviving trajectories at time t is given by 



introducing the dimensionless time r — Xah ^ t, we find 

o ma 2 ' 



that 



.5 2 . 3tt 
l-P(t) = cxp^--r+-sm T r 



(14) 



Next, we consider a Gaussian-like wave packet of free 
particles traveling toward an absorbing wall at x = 
with positive mean velocity kg. That is, in order to main- 
tain the zero boundary condition on the wall the initial 
wave function is the difference between two antisymmet- 
ric Gaussians relative to the absorbing wall. It follows 
that 



_0_ 

dx 



K(0,t] 



a h 2 



16tt 2 



exp 



a 2 (a 2 . 



16 



t 2 \2 

Ira 2 ) 



(15) 



16 



t 2 
Am 2 



3 



hence 

/ \^'(Q,t)\ 2 dt < oo. (16) 
Jo 

Thus 

R = lim [1 - P(t)] > 0, 

t — »oo 

that is, the wave packet is only partially absorbed. This 
means that the "reflected" wave consists of trajectories 
that turned around before propagating into the absorb- 
ing wall where absorption occurs. The discount of the 
wave function occurs when the packet is at the wall, as 
can be seen from eqs. (JL5|) and (|l6|). Thus R plays the 
role of a reflection coefficient. This is neither the usual 
reflection coefficient for a finite potential barrier nor that 
for an infinite barrier. 

The two examples can be combined into a simple ex- 
perimental setup of a cavity with absorbing walls and 
an absorbing detector at one end. A particle travelling 
along the axis of the cavity fits the first example in the 
the transverse direction and the second example in the 
direction of the cavity axis. Thus the decay law is the 
product of the two decay laws described above. Further 
examples and applicatins are discussed in g . 

Discussion 



complex valued and time dependent. This potential de- 
pends on the initial energies of the system [Q. Thus, 
in our formalism, the trajectories are given an actual 
physical interpretation as the possible trajectories of a 
quantum particle. This is analogous to the trajectory 
approach to diffusion in probability theory HnlO] . 

Quantum mechanics without absorption is recovered 
from our formalism when the absorbing boundaries are 
moved to infinity. In higher dimensions, quantum me- 
chanics without absorption can be recovered from our for- 
malism by putting absorbing regions with variable den- 
sity in a half space, say. As the density increases, the 
boundary of the half space becomes a totally absorbing 
wall and as the density decreases to zero, quantum me- 
chanics is recovered. 

The examples demonstrate the expected phenomenon 
that particles that reach the absorbing boundary are 
partially reflected and partially absorbed. In either case 
the decay pattern of the wave function seems to be new. 
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Absorption in a surface is different than absorption in 
the bulk across the surface in that Feynman trajectories 
do not propagate across the surface in the former but 
do in the latter case. This letter is concerned with ab- 
sorption in a surface. Absorption in the bulk requires a 
separate theory. The basic assumption in our model is 
that Feynman trajectories that propagate into the sur- 
face are instantaneously absorbed and the probability of 
the remaining trajectories is discounted by the probabil- 
ity of the absorbed trajectories at each time step. Thus 
the instantaneous discount factor is proportional to the 
probability density at the surface at each time step. The 
proportionality constant, denoted A, is a characteristic 
length, in analogy with the scattering length, the mean 
free path 0, or a typical Compton wavelength. It serves 
as a fudge parameter in this theory and is expected to 
be a measurable quantity. It may depend on the en- 
ergy of the particles, on the temperature of the absorbing 
medium, and so on. 

Our derivation does not start with a Hamiltonian, but 
rather with an action of bounded trajectories. The re- 
sulting decaying wave function corresponds to a classical 
quantum system with a Hamiltonian whose potential is 
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